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HIGHLIGHTS 


•  A  new  approach  to  the  estimation  of  PEM  fuel  cell  impedance  is  proposed. 

•  The  approach  employees  PRBS  waveform  as  perturbation  signals. 

•  The  impedance  characteristic  is  computed  based  on  continuous  wavelet  transform. 

•  Guidance  for  optimal  parameter  selection  is  presented. 

•  The  impedance  characteristic  of  a  PEM  fuel  cell  is  obtained  in  60  seconds. 


ARTICLE 


N  F  O 


A  B  S  T  R 


C  T 


Article  history: 

Received  29  August  2013 
Received  in  revised  form 
9  December  2013 
Accepted  17  December  2013 
Available  online  30  December  2013 


Keywords: 

PEM  fuel  cell 
Impedance  characteristic 
Electrochemical  impedance  spectroscopy 
Pseudo-random  binary  sequence 
Continuous  wavelet  transform 
Morlet  wavelet 


This  paper  proposes  an  approach  to  the  estimation  of  PEM  fuel  cell  impedance  by  utilizing  pseudo¬ 
random  binary  sequence  as  a  perturbation  signal  and  continuous  wavelet  transform  with  Morlet 
mother  wavelet.  With  the  approach,  the  impedance  characteristic  in  the  frequency  band  from  0.1  Hz  to 
500  Hz  is  identified  in  60  seconds,  approximately  five  times  faster  compared  to  the  conventional  single¬ 
sine  approach.  The  proposed  approach  was  experimentally  evaluated  on  a  single  PEM  fuel  cell  of  a  larger 
fuel  cell  stack.  The  quality  of  the  results  remains  at  the  same  level  compared  to  the  single-sine  approach. 

©  2013  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Electrochemical  Impedance  Spectroscopy  (E1S)  is  a  widely 
adopted  method  for  the  characterization  and  diagnostics  of  Proton 
Exchange  Membrane  (PEM)  fuel  cells.  Its  use  is  comprehensively 
covered  in  reviews  [1—4],  Nowadays,  for  the  purpose  of  the  EIS,  the 
impedance  characteristic  is  measured  at  a  fixed  operating  point 
defined  by  current  Idc-  To  perturb  the  cells,  small-amplitude  sinu¬ 
soidal  perturbations  with  frequencies  w  e  Q  are  applied  [5], 
Assuming  linearity  of  the  fuel  cell  dynamics  around  the  operating 
point  /dc,  the  fuel  cell  responds  with  sinusoid  voltage  change  with 
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the  same  frequency  w  and  particular  amplitude  and  phase.  From 
these  signals  one  can  compute  the  fuel  cell  impedance  at  the  fre¬ 
quency  to.  In  order  to  cover  the  whole  desired  frequency  range  Q 
one  has  to  repeat  the  same  procedure  multiple  of  times.  Such  an 
approach  can  be  time  consuming  and  in  some  cases  even  inap¬ 
propriate  for  on-line  monitoring  of  operational  fuel  cell  system. 
Therefore,  to  speed  up  the  experiment  with  minimal  intrusion  into 
system’s  operation,  a  new  time  efficient  approach  for  measurement 
of  impedance  characteristic  is  proposed  in  this  paper.  The  method 
uses  Pseudo-Random  Binary  Sequence  (PRBS)  as  perturbation 
signal  and  the  impedance  is  estimated  using  continuous  wavelet 
transform  (CWT)  with  Morlet  mother  wavelet. 

With  currently  adopted  single-sine  perturbation  approach,  the 
impedance  plot  of  a  fuel  cell  is  constructed  at  discrete  frequencies. 
The  approach  provides  precise  impedance  estimates  because  all  the 
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energy  of  the  perturbation  signal  is  condensed  at  one  frequency 
only,  as  such  enabling  maximal  perturbation  at  this  frequency  at 
the  expense  of  long  execution  time  period  of  the  measurements. 
Despite  the  obvious  deficiencies,  the  approach  has  been  employed 
by  the  majority  of  authors.  For  instance  Fouquet  et  al.  [6]  used  the 
EIS  based  methodology  to  monitor  the  state-of-health  of  a  PEM  fuel 
cell,  estimating  the  impedance  at  frequencies  ranging  from  0.1  Hz 
to  1  kHz  with  10  frequency  points  per  decade.  Reportedly,  one 
frequency  sweep  took  approximately  5  minutes  to  complete. 
Similarly,  Le  Canut  et  al.  [7]  successfully  detected  water  manage¬ 
ment  faults  and  catalyst  poisoning  occurrence  by  examining  the 
impedance  in  a  higher  frequency  band  from  1  Hz  to  10  kHz,  where 
one  full  measurement  sweep  took  approximately  3  minutes.  Long 
measurement  times  lead  to  the  problems  associated  with  low- 
frequency  stability  and  non-linearity  of  a  fuel  cell.  Additionally, 
from  a  practical  point  of  view,  such  time  consuming  measurements 
make  this  approach  inadequate  for  on-line  industrial  applications. 

Some  attempts  to  reduce  the  required  time  of  impedance 
measurement  have  already  been  made.  Brunetto  et  al.  [8]  proposed 
a  multi-sine  approach,  where  the  fuel  cell  under  test  is  perturbed 
with  a  signal  composed  of  multiple  sinusoids,  hence  economizing 
the  required  measurement  time.  In  the  work  of  Wasterlain  et  al.  [9], 
the  authors  clearly  exposed  the  idea  of  employing  multi-sine  sig¬ 
nals  in  their  further  work,  even  though  the  single-sine  approach 
was  used  in  the  work.  In  spite  of  the  fact  that  the  multi-sine 
approach  shortens  the  required  time  period  for  the  EIS  measure¬ 
ments,  the  approach  is  not  widely  adopted.  One  possible  reason 
might  be  that  time  reduction  is  not  significant  enough  to  justify 
additional  complexity  of  measurement  equipment  and 
computation. 

Addressing  the  issues  of  time  consuming  measurements  and 
complexity  of  measurement  equipment,  this  paper  presents  a  novel 
approach  for  fast  measurement  of  the  PEM  fuel  cell  impedance  for 
the  EIS  purposes.  The  proposed  approach  utilizes  PRBS  as  a 
perturbation  signal.  The  PRBS  is  well-established  perturbation 
signal  in  the  field  of  system  identification  since  it  has  suitable 
white-noise-like  spectral  properties  in  a  predefined  frequency 
band  [10],  However,  to  the  best  of  the  authors  knowledge  PRBS  has 
been  neglected  for  EIS  impedance  measurements  in  the  context  of 
fuel  cell  diagnostics.  Recently  published  study  by  Boghani  et  al.  [11  ] 
is  the  nearest  available  study  to  the  field  of  “conventional”  fuel  cells 
that  uses  the  PRBS  system  identification  concept.  In  their  study,  the 
PRBS  was  utilized  together  with  step  waveform  to  analyse  perfor¬ 
mance  of  a  microbial  fuel  cell  by  measuring  the  time  constant  and 
steady-state  gain  of  a  predefined  first  order  model  over  the  oper¬ 
ation  range  of  a  cell.  In  the  same  manner,  Fairweather  et  al.  [12] 
employed  a  very  low  frequency  PRBS  in  order  to  estimate  model 
parameters  of  a  valve-regulated  lead-acid  batteries  affected  by 
operational  temperature.  Unlike  these  studies,  the  novel  approach 
presented  in  this  paper  deals  with  measuring  the  impedance 
characteristic  in  a  frequency  band  from  0.1  Hz  to  500  Hz  and  it  is 
not  limited  or  predefined  by  any  linear  or  non-linear  model. 

The  main  advantage  of  the  proposed  PRBS  approach  compared 
to  the  single-sine  and  multi-sine  ones  is  a  shorter  time  period 
required  for  the  measurements.  In  the  case  of  0.1  Hz— 500  Hz  fre¬ 
quency  band,  the  measurements  take  60  seconds  to  complete.  The 
length  of  measurement  time  is  tightly  linked  to  the  lower  frequency 
limit  of  the  frequency  band,  therefore  only  slightly  higher  limit  (e.g. 
0.5  Hz  instead  of  0.1  Hz)  significantly  shortens  the  required  mea¬ 
surement  time  period.  Additionally,  as  the  cell  is  perturbed  by  the 
broadband  signal,  the  value  of  impedance  can  be  computed  at  any 
frequency  of  interest  just  with  properly  adjusted  parameters  of  the 
computational  algorithm. 

In  the  case  of  sinusoid  waveform  signals,  the  Fourier  transform 
is  applied  to  transform  the  signals  into  frequency  domain  [8,13], 


However,  when  using  broad-band  random  signals  such  as  PRBS,  it  is 
preferred  to  perform  time-frequency  analysis  for  the  impedance 
computation.  One  approach  is  the  application  of  Short-time  Fourier 
transform  or  wavelet  transform.  The  former  one  suffers  from  fixed 
time-frequency  resolution,  whereas  the  latter  one  allows  flexibility 
in  the  time— frequency  resolution  by  achieving  good  time  resolu¬ 
tion  for  high-frequency  events,  and  good  frequency  resolution  for 
low-frequency  event.  Therefore,  in  this  study  CWT  approach  based 
on  the  Morlet  mother  wavelet  was  employed  for  impedance 
computation.  With  properly  tuned  CWT  parameters,  this  approach 
provides  reliable  impedance  results  along  the  entire  frequency 
band.  Additionally,  it  yields  statistical  information  about  confi¬ 
dence  interval  of  the  impedance  measurement. 

The  wavelet  transform  is  widely  used  tool  in  a  variety  of  sci¬ 
entific  fields  ranging  from  diagnostics  of  human  cardiovascular 
systems  [14]  and  analysis  of  life  signals  in  mammals  [15]  to  fault 
detection  in  mechanical  drives  [16],  In  the  field  of  PEM  fuel  cell 
diagnostics,  Steiner  et  al.  [17]  used  discrete  wavelet  transform  to 
diagnose  PEM  fuel  cell  flooding  directly  from  the  voltage  signals. 
The  flooding  was  diagnosed  based  on  the  pattern  of  the  wavelet 
packet  coefficients  derived  with  the  wavelet  packet  transform. 

The  paper  also  addresses  the  problem  of  proper  parameter  se¬ 
lection  of  both  PRBS  and  CWT  with  the  Morlet  wavelet,  which  has 
significant  influence  on  the  effectiveness  of  the  proposed  approach. 
Special  attention  is  given  to  the  computationally  efficient  imple¬ 
mentation  of  the  algorithms  which  brings  significant  practical  ad¬ 
vantages  in  implementation  of  on-line  EIS  based  diagnostic 
systems.  Finally,  the  paper  presents  experimental  results,  where 
the  PRBS  approach  was  applied  to  a  commercial  PEM  fuel  cell 
system.  The  PEM  fuel  cell  impedance  results  acquired  using  the 
PRBS  approach  are  thoroughly  discussed.  The  impedance  charac¬ 
teristic  measured  with  use  of  PRBS  perturbation  signal  is  compared 
to  the  reference  impedance  characteristic  obtained  by  standard 
single-sine  approach. 


2.  Pseudo-Random  Binary  Sequence 


According  to  Ljung  [10],  PRBS  is  a  periodic  and  deterministic 
signal  with  white-noise-like  properties.  As  such,  it  is  highly 
appropriate  as  an  input  test  signal  for  system  identification  pur¬ 
poses.  The  maximum  length  PRBS  can  be  generated  by  employing 
feedback  shift  registers  [18,19],  The  maximum  length  PRBS  is 
characterized  by  its  order  n,  the  maximum  length  period  N, 
amplitude  a  and  its  sampling  period  At.  The  order  n  determines  the 
maximum  length  period  N  (number  of  discrete  points  in  time, 
when  the  PRBS  signal  can  change  its  value): 


The  discrete  power  spectrum  density  4>  of  maximum  length 
PRBS  signal  with  amplitude  a  and  length  N  is  [19]: 


ri  f  f  m  —  0 

bd(m)  =  < 

\ 0  <m<N 


Therefore,  the  power  spectrum  of  periodical  PRBS  with  period  N 
and  sampling  period  At  is: 


(3) 


Fig.  1  resembles  the  power  spectrum  of  the  PRBS  signal  given  by 
the  Equation  (3).  It  can  be  seen,  that  the  value  of  the  power  spec¬ 
trum  hits  zero  exactly  at  integer  multiples  of  sampling  frequency /s, 
which  is  related  to  the  sampling  period  At  as  follows: 
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i p(t)  =  TC-1/4  (e,6,ot  -  e-"o/2)e-ta/2,  (8) 

where  wo  stands  for  the  ratio  between  the  highest  and  the  second 
most  highest  peak  and  is  usually  set  to  wo  =  5.  For  such  wo  >  5  the 
second  term  in  (8)  can  be  neglected. 

The  Morlet  wavelet  (8)  is  an  analytical  function  i.e.  it  has  only 
positive  frequencies.  Consequently,  the  wavelet  coefficients  Wf[s,u ) 
in  (6)  are  complex  values.  Therefore,  the  CWT  is  plotted  in 
modulus-phase  form  i.e.  one  plot  for  |  (f.i/'u.s)  |  and  another  plot  for 
argCtfi/ys)).  both  on  time-frequency  plane.  As  a  result,  at  each  time 
translation  u  and  scale  s  the  CWT  with  Morlet  wavelet  gives  the 
instantaneous  amplitude  and  phase. 

3.2.  Scale  to  frequency  transform  for  Morlet  wavelet 


Fig.  1.  Power  spectrum  of  the  PRBS  signal. 


fs 


At 


(4) 


However,  the  system  under  test  has  to  be  sufficiently  perturbed, 
therefore  the  useful  frequency  band  is  determined  by  -3  dB  fre¬ 
quency  limit.  In  this  way,  the  useful  frequency  band  /b  is  approxi¬ 
mately  determined  as: 


The  frequencies  in  CWT  (6)  are  represented  by  scale  s,  which 
corresponds  to  a  particular  dilatation  of  the  mother  wavelet 
\p{x).  Transforming  scale  to  frequency  for  the  Morlet  wavelet  is 
defined  as  [22]: 


O) 


where  s  is  the  scale  of  CWT  and  wo  is  the  parameter  of  the  Morlet 
wavelet  governing  the  peak  ratios. 


h  =  ys 


(5)  3.3.  Cone  of  influence 


3.  Continuous  wavelet  transform 

The  application  of  Fourier  transform  is  suitable  for  roughly 
stated  smooth  periodic  signals.  However  for  signals  describing 
transient  phenomena,  located  in  a  narrow  time  interval,  the 
application  of  Fourier  transform  turns  to  be  inadequate.  Moreover, 
the  Fourier  transform  performs  poorly  at  lower  frequencies. 

Wavelet  transform  is  based  on  a  set  of  specifically  designed 
functions  called  wavelets.  The  continuous  wavelet  transform  (CWT) 
of  a  square  integrable  function /(t)eL2(R)  is  defined  as  [20]: 

Wf(s,u)  =  <f(t),t u,s(t)>  =  J  /((#*, s(t)dt,  (6) 

where  i//u,s (t)  is  a  scaled  and  translated  version  of  the  mother 
wavelet  i f{t): 


The  CWT  transform  (6)  describes  the  analysed  signal  J[t)  on  the 
time-scale  plane.  The  conversion  between  scale  s  and  actual  fre¬ 
quency  /  is  straightforward  and  depends  on  the  selection  of  the 
mother  wavelet  (7). 

Each  scaling  s  alters  the  time-frequency  resolution  of  the 
mother  wavelet.  As  a  result,  CWT  offers  adaptive  time-frequency 
resolution.  High  frequency  resolution  is  preserved  for  lower  fre¬ 
quencies,  whereas  high  time  resolution  is  preserved  for  higher 
frequencies. 

3.1.  The  Morlet  wavelet 

For  the  impedance  estimation,  one  needs  information  about  the 
instantaneous  amplitude  and  phase  of  the  electrical  current  i(t)  and 
voltage  u(t).  Therefore,  a  logical  choice  for  mother  wavelet  (7)  is  the 
Morlet  wavelet  defined  as  [21]: 


Due  to  the  finite  length  of  the  measured  signals,  the  wavelet 
coefficients  (6)  will  be  influenced  by  the  discontinuities  at  the 
beginning  and  the  end  of  the  signal.  These  discontinuities  intro¬ 
duce  false  broad  frequency  components.  The  presence  of  these  ef¬ 
fects  is  different  over  different  wavelet  scales,  defining  the  so-called 
cone  of  influence.  It  is  defined  as  the  e-folding  time  for  the  auto¬ 
correlation  of  wavelet  power  at  each  scale  [22],  which  for  the 
Morlet  wavelet  reads: 

e  =  sV 2.  (10) 

The  cone  of  influence  for  the  Morlet  wavelet  is  shown  as 
hatched  region  in  Fig.  2.  The  hatched  regions  show  the  wavelet 
coefficients  whose  values  are  influenced  by  the  edge  effects.  These 
coefficients  should  be  omitted  from  further  analysis.  In  Fig.  2,  it  can 
be  noticed  that  all  of  the  low  frequency  wavelet  coefficients  are 
influenced  by  the  discontinuities  occurring  on  both  edges.  One  can 
overcome  this  limitation  by  either  using  longer  signal  or  by  over- 
sampling  the  existing  one. 

3.4.  Influence  of  the  Morlet  wavelet  central  frequency 

The  definition  of  Morlet  wavelet  (8)  states  that  the  wavelet  is 
dependent  upon  parameter  oiq.  The  choice  of  this  parameter  is  a 
compromise  between  localization  in  time  and  frequency  [14,15], 
For  smaller  wo,  the  shape  of  the  wavelet  favours  localization  of 
singular  time  events,  whilst  for  larger  w0  more  periods  of  the  sine- 
carrier  in  the  Gaussian  window  (envelope)  make  the  frequency 
localization  better. 

The  broad  frequency  band  and  relatively  short  duration  of 
measured  signals  (in  case  of  lower  frequencies)  raise  a  problem  in 
relation  to  variance  of  the  calculated  impedance.  In  order  to 
minimize  the  variance  the  frequency  resolution  is  increased  by 
increasing  the  central  frequency  of  the  wavelet  wo. 

While  this  minimizes  the  variance  for  high  frequencies  it 
simultaneously  raises  the  variance  for  lower  frequencies.  This  is 
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Fig.  2.  Cone  of  influence  for  the  Morlet  wavelet  with  edge  e  defined  by  Equation  (10). 


due  to  the  fact  that  at  lower  frequencies  the  whole  signal  contains 
less  periods  than  the  wavelet  itself,  which  increases  the  variance.  To 
remedy  this  situation  the  parameter  w0  is  lowered  to  decrease  the 
number  of  periods  in  the  wavelet.  While  this  trade-off  lowers  fre¬ 
quency  resolution  it  improves  the  reliability  of  the  calculated 
impedance. 


4.  Experimental 

4.1.  Experimental  setup 

The  experimental  setup  consisted  of  a  PEM  fuel  cell  power  unit, 
electronic  load,  arbitrary  function  generator,  voltage  and  current 
measurement  equipment,  and  desktop  computer.  The  measure¬ 
ments  were  performed  on  PEM  fuel  cells  of  a  commercial  8.5  kW 
fuel  cell  system  HyPM  HD8-200,  produced  by  Hydrogenics  Inc.  The 
stack  of  the  system  consists  of  80  PEM  fuel  cells  with  surface  area  of 
200  cm2  each. 

The  electronic  load  was  connected  to  the  fuel  cell  system  and  it 
was  used  to  perturb  the  system  with  PRBS  perturbation  signal  in 
galvanostatic  mode.  The  waveform  of  the  load  current  was 
controlled  by  an  arbitrary  function  generator.  The  generator  was 
generating  the  PRBS  reference  signal,  which  was  feed  to  the  elec¬ 
tronic  load  as  a  reference  signal  for  the  load  current.  Block  diagram 
in  Fig.  3  presents  the  experimental  setup.  Here,  all  the  signal  flows 
can  be  distinctly  observed. 
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Fig.  3.  Block  diagram  of  the  experimental  setup. 


Voltage  and  current  signals  were  measured  with  an  in-house 
developed  fuel  cell  voltage  monitor  (FCVM)  described  in  Ref.  [23] 
and  presented  in  Fig.  4.  The  top  view  of  the  figure  shows  the 
printed  circuit  board  of  the  FCVM.  The  board  is  mounted  directly  on 
top  of  the  stack  and  the  galvanic  contacts  to  all  the  fuel  cells  are 
implemented  by  spring  contact  probes,  which  are  visible  in  the  side 
view  of  Fig.  4.  The  FCVM  provides  means  for  fast  and  accurate 
measurements  of  voltage  and  current  signals  of  any  individual  fuel 
cell  inside  the  stack  (the  contacts  to  all  fuel  cells  of  the  stack  are 
clearly  seen  at  the  bottom  edge  of  Fig.  4).  The  FCVM  utilizes  a  Hall- 
effect  current  sensor  to  measure  current,  and  an  analogue  circuitry 
to  precisely  measure  changes  in  voltage  of  any  individual  fuel  cell 
inside  the  stack.  Any  cell  of  the  stack  can  be  selected  to  be 
measured,  however  only  three  of  them  can  be  measured  simulta¬ 
neous  at  the  same  time. 

With  the  FCVM,  voltage  and  current  signals  are  measured  with 
the  resolution  of  80  pV  and  10  mA,  respectively,  at  sampling  fre¬ 
quency  of  5  kHz.  The  FCVM  has  the  frequency  range  up  to  660  Hz. 
The  acquired  measurements  are  transferred  to  a  PC  via  CAN  bus 
interface.  The  signal  processing  and  impedance  computational  al¬ 
gorithms  are  performed  on  the  PC. 

4.2.  Selection  of  the  PRBS  bandwidth  and  Morlet  wavelet  central 
frequency  values 

Due  to  the  limitations  set  by  the  FCVM  (i.e.  frequency  range  of 
660  Hz),  the  impedance  characteristic  was  measured  in  the  fre¬ 
quency  band  from  0.1  Hz  to  500  Hz.  Since  such  a  frequency  band 
spans  over  more  than  three  frequency  decades,  two  PRBS  wave¬ 
forms  with  different  sampling  frequency  fs  were  employed  as 
perturbation  signals.  First  one  with  fs  —  500  Hz  covered  the  fre¬ 
quencies  from  0.1  Hz  to  10  Hz,  whereas  the  other  with  /s  =  2  kHz 


Fig.  4.  Top  and  side  view  of  the 
PEM  fuel  cells. 


FCVM  connected  i 
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covered  the  frequencies  from  10  Hz  to  500  Hz.  According  to  (5),  the 
useful  bandwidths  ft  of  the  two  signals  are  approximately  166  Hz 
and  660  Hz,  respectively. 

Accordingly  to  the  analysis  given  in  Section  3.4,  for  proper 
estimation  of  impedance  at  various  frequencies,  the  central  fre¬ 
quency  wq  of  the  used  Morlet  wavelet  was  accordingly  adjusted 
along  the  frequency  axis.  Three  values  of  the  central  frequency  wo 
were  selected,  each  covering  separate  frequency  band. 

Fig.  5  graphically  depicts  the  selected  bandwidth /b  of  the  PRBS 
signal  and  the  corresponding  central  frequency  wo  of  the  Morlet 
wavelet  for  different  frequency  ranges.  The  selection  of  the  PRBS 
bandwidth  was  done  accordingly  to  the  analysis  given  in  Section  2. 
The  wavelet  central  frequency  wo  was  chosen  in  line  with  the 
analysis  given  in  Section  3.4. 

4.3.  Experimental  procedure 

During  the  experiments,  the  fuel  cell  system  was  operating  at 
constant  operating  and  environmental  conditions.  The  DC  current 
operating  point  ftc  was  set  to  50  A,  whereas  the  amplitude  of  the 
superposed  PRBS  waveform  was  set  to  2  A.  As  such,  the  peak-to- 
peak  amplitude  was  4  A  and  therefore  small  enough  not  to  cause 
difficulties  due  to  non-linearity  of  the  PEM  fuel  cells. 

The  order  n  of  the  PRBS  was  selected  to  be  n  =  16.  Therefore, 
according  to  ( 1 )  the  length  of  PRBS  sequence  is  65535  samples.  The 
two  PRBS  signals  were  applied  to  the  fuel  cell  system  in  consecutive 
manner  one  after  another.  The  first  PRBS  signal  (ft  =  166  Hz)  was 
applied  for  45  s  and  the  second  (ft  =  660  Hz)  for  15  s,  therefore  one 
individual  measurement  of  impedance  characteristic  took  60  s  in 
total  to  complete.  The  first  PRBS  signal  was  applied  for  longer  time 
period  since  it  excites  lower  frequencies.  Fig.  6  shows  time  plots  of 
current  and  voltage  signals  of  a  PEM  fuel  cell  during  measurements. 
Current  was  measured  as  an  absolute  value,  whereas  the  fuel  cell 
voltage  was  measured  as  differential  value. 

5.  Results  and  discussion 

Utilizing  the  PRBS  perturbation  signals  and  CWT  computational 
algorithm,  the  obtained  impedance  characteristic  of  a  PEM  fuel  cell  is 
presented  in  Fig.  7  in  the  form  of  Bode  plot.  The  impedance  char¬ 
acteristic  was  computed  at  450  different  frequency  points.  However, 
since  the  measured  current  and  voltage  signals  comprise  all  fre¬ 
quency  components  in  the  frequency  band,  the  algorithm  allows  to 
arbitrary  choose  frequencies  and  the  number  of  frequency  points  at 
which  the  impedance  is  to  be  computed.  In  every  such  a  configu¬ 
ration  set  of  frequency  points,  the  input  signals  for  the  algorithm  are 
the  same,  the  only  thing  is  the  computation  time  (i.e.  more  frequency 
points  lead  to  longer  computation  time,  and  vice  versa). 

The  plotted  impedance  characteristic  in  Fig.  7  is  the  median 
value  of  the  impedance  through  measurement  time  period. 
Alongside  the  median  impedance  value,  the  10th  and  90th  per¬ 
centiles  are  provided  in  order  to  show  the  reliability  of  the 
impedance  estimation.  As  such,  the  two  percentiles  account  for  80% 
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Fig.  5.  Sampling  frequency  fs  of  applied  PRBS  signals  and  utilized  central  frequencies 
w0  of  complex  Morlet  wavelet  in  relation  to  frequency  scale  /  =  [0.1,  500]  Hz. 


Fig.  6.  Measured  current  and  voltage  signals  during  the  experiment. 


of  the  population.  Additionally,  the  figure  also  provides  the  refer¬ 
ence  values  of  the  impedance,  which  are  presented  with  dots.  The 
reference  values  were  measured  with  single-sine  technique  at  4 
points  per  decade. 

Fig.  7a  shows  the  Bode  amplitude  plot  of  the  measured  PEM  fuel 
cell  impedance.  In  the  frequency  band  from  1  Hz  to  500  Hz,  the 
results  exhibit  high  alignment  with  the  reference  impedance 
values.  Moreover,  the  span  between  10th  and  90th  percentile  is 
quite  narrow,  meaning  that  the  impedance  results  exhibit  high 
confidence  in  this  frequency  band.  At  slightly  lower  frequencies  in 
the  frequency  band  between  0.5  and  1  Hz,  the  median  amplitude 
value  of  the  impedance  still  exhibit  steady  behaviour,  but  the  80% 
confidence  interval  is  noticeable  extended.  In  the  lowest  frequency 
band  from  0.1  to  0.5  Hz,  the  80%  confidence  interval  is  wide, 
therefore  the  results  in  this  frequency  band  has  always  to  be 
interpreted  together  with  the  confidence  interval.  However,  the 
median  value  is  informative,  but  it  should  not  be  declared  as  the 
exact  value  of  the  impedance  amplitude  without  considering  the 
confidence  interval.  The  reasons  for  the  occurrence  of  result  with 
such  relatively  high  uncertainty  can  be  found  in  at  least  two  causes. 
Firstly,  impedance  estimation  at  low  frequencies  requires  suffi¬ 
ciently  long  signals  in  order  to  obtain  more  reliable  results.  Sec¬ 
ondly,  the  aforementioned  cone  of  influence  of  the  Morlet  wavelet 
decreases  the  number  of  eligible  samples  for  calculating  the 
impedance  statistics  at  very  low  frequencies. 

The  Bode  phase  plot  of  the  measured  PEM  fuel  cell  impedance  is 
presented  in  Fig.  7b.  The  previous  discussion  about  amplitude  re¬ 
sults  applies  equally  to  the  phase  results.  As  can  be  seen  in  Fig.  7b, 
the  results  indicate  that  the  measured  PEM  fuel  cell  exhibits 
inductive  behaviour  at  low  frequencies.  According  to  Roy  and 
Orazem  [24],  the  low-frequency  inductive  behaviour  can  be 
attributed  to  PEM  fuel  cell  characteristic  itself  and  not  solely  to  non¬ 
stationary  artefacts.  Roy  et  al.  [25]  suggested  the  reactions  of 
hydrogen  peroxide  intermediate  formation  and  PtO  formation  as 
possible  causes  of  the  low-frequency  inductive  behaviour.  Since  the 
inductive  behaviour  was  observed  with  the  PRBS  approach  as  well 
as  with  the  single-sine  one,  most  probably,  the  inductive  behaviour 
was  not  introduced  by  the  PRBS  or  CWT  methodology. 

Fig.  8  summarizes  the  impedance  result  in  the  form  of  Nyquist 
plot.  The  plot  presents  only  the  median  value  of  the  impedance 
without  80%  confidence  interval  since  there  is  no  straightforward 
manner  of  presenting  the  confidence  interval  of  the  impedance 
characteristic  in  the  Nyquist  plot. 

The  Nyquist  plot  in  Fig.  8  clearly  expresses  two  semicircle  arcs  in 
the  PEM  fuel  cell  impedance  characteristic.  The  high-frequency  arc 
ranges  from  approximately  5  Hz  up  to  high  frequencies.  The  low- 
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Bode  phase  plot. 


frequency  one  starts  at  approximately  0.5  Hz  and  extends  up  to 
5  Hz,  where  the  high-frequency  arc  starts.  The  Nyquist  plot  also 
clearly  shows  the  already  discussed  inductive  behaviour  at  lowest 
frequencies. 

The  study  has  confirmed  the  hypothesis  that  the  measurements 
of  the  PEM  fuel  cell  impedance  utilizing  the  PRBS  perturbation 
signal  can  be  performed  in  a  substantially  shorter  time  period  than 
with  accepted  single-sine  approach.  In  the  case  of  the  proposed 


Re(Z)  [mSl] 


Fig.  8.  Nyquist  plot  of  the  measured  PEM  fuel  cell  impedance. 


PRBS  approach,  the  required  time  period  for  the  measurements  is 
approximately  five  times  shorter  compared  to  the  single-sine 
technique.  With  such  a  reduction  of  time  required  for  the  mea¬ 
surements,  the  feasibility  and  acceptability  of  the  EIS  based  di¬ 
agnostics  in  industrial  applications  is  significantly  enhanced. 
Moreover,  unlike  the  single-sine  approach,  where  the  impedance 
can  be  obtained  only  at  frequencies  that  have  been  perturbed,  in 
the  case  of  the  proposed  PRBS  approach,  the  impedance  can  be 
computed  at  any  frequency  of  interest  by  simple  adjustment  of  the 
computational  algorithm. 

6.  Conclusion 

The  PEM  fuel  cell  impedance  characteristic  obtained  by 
employing  PRBS  perturbation  signals  and  CWT  computation  algo¬ 
rithm  shows  high  similarity  to  the  impedance  characteristics  ob¬ 
tained  by  the  single-sine  approach.  The  similarity  between  the  two 
characteristics  confirm  the  suitability  of  the  proposed  PRBS 
perturbation  and  CWT  computational  approach.  Regardless  of  the 
fact,  that  the  present  study  has  been  performed  on  a  PEM  fuel  cell, 
the  proposed  approach  for  impedance  measurement  can  be  easily 
applied  to  other  types  of  electrochemical  cells  as  well. 

In  this  study  impedance  characteristic  in  the  frequency  band 
from  0.1  Hz  to  500  Hz  was  measured.  For  such  a  measurement,  the 
PRBS  approach  takes  60  s  to  complete.  However,  if  the  impedance 
only  at  frequencies  above  1  Hz  is  to  be  measured,  the  required  time 
for  the  measurements  is  further  reduced.  Additionally,  the  upper 
bandwidth  frequency  can  be  freely  selected  even  at  higher 
frequencies. 

The  computational  algorithm  besides  the  median  value  of  the 
population  provides  also  the  confidence  interval  of  the  results.  The 
information  about  the  confidence  interval  is  very  beneficial,  since  it 
expresses  how  far  the  results  can  be  trusted.  In  this  manner,  the 
study  revealed  that  the  impedance  results  below  0.5  Hz  expresses 
wide  confidence  interval,  therefore  they  have  to  be  interpreted 
with  great  care  if  used  for  the  EIS  based  diagnostics  or  are  employed 
for  any  other  purposes. 

The  presented  approach  of  measuring  PEM  fuel  cell  impedance 
clearly  provides  the  impedance  characteristic.  However,  the 
approach  still  uses  some  kind  of  perturbation  signal.  Therefore,  the 
possibility  of  using  no  external  perturbation  signal  at  all  should  be 
investigated  in  the  future  work.  This  way,  solely  the  readily  avail¬ 
able  signals  (i.e.  the  operation  signal  of  the  fuel  cell  system)  should 
be  employed  for  impedance  characteristic  identification.  Such  ob¬ 
tained  data  could  be  further  used  for  improving  the  operation  of  in- 
the-field  fuel  cell  power  units  [26], 

In  this  study,  the  approach  was  applied  to  estimate  the  imped¬ 
ance  of  a  PEM  fuel  cell.  However,  the  employed  hardware  and 
proposed  approach  are  broadly  applicable  in  the  field  of  electro¬ 
chemical  impedance  spectroscopy.  With  minor  modifications  the 
approach  may  be  applied  to  any  other  type  of  electrochemical  cells, 
for  instance  lithium  ion  cells. 
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